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Abstract 

We construct PT-symmetric quantum mechanical models with an 0(iV)-symmetric interaction 
term of the form —g (x 2 ) 2 /N. Using functional integral methods, we find the equivalent Hermitian 
model, which has several unusual features. The effective potential obtained in the large- TV limit of 
the Hermitian form is shown to be identical to the form obtained from the original PT-symmetric 
model using familiar constraint field methods. The analogous constraint field prescription in four 
dimensions suggests that —g ( 4> 2 \ /N PT-symmetric scalar field theories are asymptotically free. 
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I. INTRODUCTION 

Since the initial discovery of PT symmetry [lj], there has been considerable progress in 
expanding both the number of PT-symmetric models and our knowledge of their proper- 
ties {2], [J Although the construction of PT-symmetric matrix models has proved to be 
relatively straightforward Q|, the construction of PT-symmetric models with fields trans- 
forming as vectors under O(N) appears to have more difficult technical issues. Nevertheless, 
the development of scalar field theory models with vector symmetry is crucial to the possible 
relevance of PT-symmetry in particle physics. On the other hand, there has been substan- 
tial progress recently in understanding the relation of the one-component —Ax 4 model to 
its equivalent Hermitian form jf], Q|, as well as work on the relation of 0(iV)-symmetric 
Hermitian models to one-component PT-symmetric models {^fl. In this paper, we extend 
this work to a PT-symmetric model with O(N) symmetry, deriving the Hermitian form for 
all values of N. We construct the large- N limit, and show that this limit can be obtained 
using the familiar method of constraint fields, but in a much less rigorous way. 

Consider a model with Euclidean Lagrangian given by 
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where g and A are non-negative. When g = 0, we have N decoupled one-dimensional systems; 
for A = 0, we have a model with O(N) symmetry. When both g and A are non-zero, the 
model has only an SV permutation symmetry. From the standpoint of PT symmetry, the 
interaction terms can be considered as members of a family of PT-invariant interactions 

N / N \ q 

j=l \ 3=1 J 

which are invariant under PT symmetry. This class of models is well-defined for p = q = 1, 
and must be defined for p, q > 1 by an appropriate analytic continuation of the Xj as 
necessary 

It is convenient to consider this model as a subset of a larger class of models, with a 
Lagrangian of the form 
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2 ( dtXj ^ 2 + \ m2x2 j ~ x ) k ^ x l ( 3 ) 

j=l 1 " j,k=l 



The classical stability of the potential for large Xj is governed by the eigenvalues of A. For 
the model of particular interest to us, 



A = XI + gP 



(4) 



where P is the projector 
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v ■■ -J 

satisfying P 2 = P. The decomposition A = A (J — P) + (g + A) P shows that A has one 
eigenvalue g + A and N — 1 eigenvalues with value A. The eigenvalue g + A is associated with 
variations in x 2 , i.e., variations in the radial direction. 



II. EQUIVALENCE TO HERMITIAN MODELS 



We will analyze the case where all eigenvalues of A are positive using functional integra 
tion. With the substitution 

x 



v 3 / Cj + iijjj 

familiar from the one-component case, becomes 



E 



— 2m 2 (cj + iijjj) 



2(cj + iipj) 

The generating function for the model is given by 
rr [#?] 
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- / dt 
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determinant and additional 



where the change of variables has generated both a functiona 
term, formally of order H 2 , in the action. As pointed out in [6j, both terms are required to 
obtain correct results in the functional integral formalism. 
The functional determinant may be written as 
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Y\ [dhj] exp I — I dt 
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which introduces a new set of fields hj. The derivation proceeds as in the single- variable case. 
After integration by parts on the hjipj terms, and adding and subtracting total derivatives, 



the functional integral over the tpj fields can be carried out exactly. The integral is both local 
and quadratic, and requires that the matrix A have positive eigenvalues for convergence. The 
result of this integration is 



Z= /n^exp/- Idt \l~J2^ 

J n I J L jk 



-h 2 + L - 2m 2 
2 



-M + h k - 2m 2 
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After discarding total derivatives, we obtain 
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which gives the Hermitian form for our general PT-symmetric model with N fields. 
In the particular case we are interested in, we have 

A- 1 = i (J - P) + -j-p. 

A g + A 



(12) 



The Lagrangian may be written as 
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(13) 



It is helpful to immediately rescale all the fields as hj — >■ v^32A/i., 
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At this point, the Sn permutation symmetry is still manifest, and it clear that the field 
YE hj plays a special role. 

In order to understand the strategy for rewriting the model in a form in which the limit 
A — > can easily be taken, it is useful to work out explicitly the case of N = 2 first. It is 
apparent that a rotation of the fields will be desirable. We define suggestively new fields a 
and 7i given by 

After some algebra and the rescaling 



(15) 



g + X 



A 



(16) 
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we arrive at 

L E = -& 2 + -7T 2 - mV - ^-n 2 + 2(g + X)<7 4 + -^-vr 4 + (8g + 12A) - 2^/gTXa. 
2 2 # + A # + A v 

(17) 

Notice the natural hierarchy between the masses for X ^ g. The 0(2) symmetric limit of 
the original PT-symmetric model is obtained in the limit A — > 0, where we have 

L E = \^y 2 + ^7T 2 - mV + 2ga 4 + 8#(j 2 7r 2 - 2y/ga. (18) 

The field 7r has no mass term, indicating its relation to the angular degrees of freedom in 
the original Lagrangian. However, radiative corrections generate a mass for the n field via 
the the cx 2 7r 2 interaction. As in the one-component case, there is a linear anomaly term, but 
only for a. 

We now turn to the more difficult case of the A — > limit for arbitrary N. As before, we 
introduce a field a defined by 

a = -L V" hj (19) 

v 3 

as well as a set of N — 1 fields Tik with k — 1, .., N — 1 related to the hj fields by a rotation 
so that a 2 + n 2 = h 2 . Each field hj can be written as 

hj = —L(7 + hi (20) 

3 3 



where ^2jhj = 0. This property is crucial in eliminating a term in L E which diverges as 
A -1 ' 2 as A — > 0. The Lagrangian now can be written as 

£ E = \^Y. \*h^ (° 2 + * 2 )+« E - (- 2 + * 2 - ^wf-^'- 

j j 

(21) 

The rescaling a — > ^(g + X) /Xa plus some careful algebra yields the A — > limit as 

= ^ + \& ~ m 2 a 2 + ^a 4 + ^Vvr 2 - ^Na (22) 

which agrees with our previous result for N = 2 , and agrees with the known result for 
a single degree of freedom if we take N — 1 and drop the n field altogether. This is a 
Hermitian form of the PT-symmetric anharmonic oscillator with O(N) symmetry, derived 
as the limit of a PT-symmetric model with Sn symmetry. The Hermitian form has several 
novel features. Note that both the and O(N) symmetries are no longer manifest, but 



there is an explicit 0(N — 1) symmetry associated with rotations of the ff field. As in the 
N = 2 case, there is no mass term for the 7? field. Furthermore, there is no (vr 2 ) 2 term, 
although there is a 7T 2 <7 2 interaction. The anomaly term again involves only ex, and breaks 
the symmetry a — > —a possessed by the rest of the Lagrangian. Analyzing the Lagrangian 
at the classical level, we see that if m 2 > , the a field is moving in a double-well potential, 
perturbed by the anomaly so that (a) > 0. On the other hand, if m 2 < 0, o moves in a 
single-well anharmonic oscillator, again with the linear anomaly term making (cr) > 0. In 
either case, the n 2 a 2 interaction will generate a mass for the n field. All of this is consistent 
with the association of a and n with the radial and angular degrees of freedom, respectively, 
in the original PT-symmetric model. 



III. LARGE- N LIMIT 



We will defer a more detailed discussion of this model for finite N, and turn to its large- N 
limit. One more rescaling a — > y/Na gives the Lagrangian 

L E = — a 2 + Ut 2 - Nm 2 a 2 + AgNa* + 16ga 2 rr 2 - N^/2g~a. (23) 

We see that the anomaly term survives in the large-iV limit, unlike the matrix model case 
J|. After integrating over the N — 1 components of the ff field, we have the large- N effective 
potential V e ff for a: 

V eff /N = -m 2 a 2 + Aga A + ^32^ 2 - y^2~ga. (24) 

It is striking that the anomaly term has virtually the same form as the zero-point energy of 
the 7? field. The anomaly term breaks the discrete a — > —a symmetry of the other terms 
of the Lagrangian, and always favors a > 0. The effective potential has a global minimum 
with a positive for m 2 > 3 2 1 / 3 g 2//3 . For m 2 < 3 2 1 / 3 g 2 / 3 , there does not appear to be a 
stable solution with a > 0, and a = is the stable solution to leading order in the 1/N 
expansion. This change in the behavior of the effective potential as m 2 is varied is not seen 
in the corresponding Hermitian model [8|, and indicates a need for care in analyzing the 
model. Based on our preliminary analysis of the Hermitian form for finite N, we believe 
that this behavior is associated with the large- N limit, and does not indicate a fundamental 
restriction on m 2 . 
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The large-iV effective potential was derived from a Lagrangian with unusual properties, 
associated with the Hermitian form of the original model. It is therefore surprising that, 
once the form of the large- N effective potential is known, it can be derived heuristically in 
a more conventional way. We start from the 0(iV)-symmetric Lagrangian 



]\[ / AT \ 2 



- (d t Xjf + -\n 2 x) 
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(25) 



and add a quadratic term in a constraint field p 



N 



2 



q I 2Np 9 Nm 2 \ . . 

L '^\r^—nr\ (26) 



yielding 



N 



4Np 2 Nm 2 p Nm 4 , x 

+ — + 27 

9 9 16^ 



If we integrate over Xj in a completely conventional way, we obtain the large- iV effective 
potential 

4p 2 m 2 p I — m 4 



9 9 160 
This is essentially identical to our previous expression after identifying p = go 2 . However, 

we lack a fundmental justification for this approach. We know that great care must be taken 

in specifying the contour of integration in typical PT-symmetric models, yet the Xj fields 

were integrated over quite conventionally. If this approach has validity, it seems likely that 

the choice of integration contours for p and x is crucial. However, only the saddle point 

matters to leading order in 1/N, so it is possible for this heuristic derivation to be correct 

even though we lack a direct, complete treatment of the original PT-symmetric model. 

IV. PT-SYMMETRIC FIELD THEORY 



If we boldly apply the constraint field approach to a PT-symmetric field theory with a 
'9 ^0 2 j interaction in d dimensions, we obtain the effective potential 

4p 2 m 2 p m 4 If d d k , r , 9 n , . 

Veff/N = — + + - / -r In \k + 8p] . 29 

ff/ 9 9 160 2 J ( 27 r) d [ } 



Models of this type were rejected decades ago [8] because of stability concerns at both 
the classical and quantum levels, although there were early indications that such theories 



were in fact sensible [9]. Within the framework of PT-symmetric models, such stability 
issues cannot be addressed without understanding the boundary conditions on the functional 
integrals. However, it is straightforward to check that renormalization of g in d = 4 gives 
an asymptotically free theory, with beta function f3 = —g 2 /2n 2 in the large-iV limit. If 
PT-symmetric scalar field theories exist in four dimensions and are indeed asymptotically 
free, the possible implications for particle physics are large, and provide ample justification 
for further work. 

Finally, we note that the construction we have used for PT-symmetric models with 0(N) 
symmetry is not the only one possible. For example, we can consider our original model 
with g > but A < 0, so that only the 0(N) symmetric term is unconventional. It would 
be interesting to know if the A — > limit of this class of models can be used to define a 
PT-symmetric model with 0(N) symmetry, and if so, if it is equivalent to the one studied 
here. 
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